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Abstract 

We evaluate external jy-emission, internal jy-emission and P^-exchange 
contributions to nonleptonic B — > D^*> decays based on the perturbative 
QCD formalism including Sudakov effects, whose ratio is found to be 1 : 
+0.2 : 0.03z at the amplitude level. We observe that the internal jy-emission 
contribution is additive to the external IV-emission contribution, and that 
the H^-exchange contribution is negligible and mainly imaginary, its real part 
being at least one order of magnitude smaller than the imaginary part. Our 
predictions are consistent with the CLEO data and with those obtained by 
the Bauer-Stech-Wirbel method. 
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1. Introduction 

The conventional approach to nonleptonic B meson decays is the Bauer- 
Stech-Wirbel (BSW) method |J, in which two parameters a x and a 2 are 
associated with the external and internal ly-emission amplitudes, respec- 
tively, and are determined by experimental data. a\ and a 2 were originally 
the linear combination of the Wilson coefficients C\ and c 2 in the effective 
Hamiltonian 

#eff = ^V cb V: d [ Cl (fi)(du)(3) + c 2 (fx)(cu)(db)} , (f) 

written as a\ = c\ + c 2 /N c and a 2 = c 2 + Ci/N c , N c being the number of 
colors. H e g was derived from the Hamiltonian 

H = ^V cb V: d (du)(cb) , (2) 

with hard gluon corrections taken into account by renormalization group 
(RG) methods. Here {q%qj) = qff^l — 75)^ represents the V — A current, Gp 
is the Fermi coupling constant, and V c b and V u d are the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements. 

The hadronic form factors involved in nonleptonic B meson decays are 
usually assumed to take, say, a monopole or dipole ansatz ||, and thus the 
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extraction of ai and 02 is model dependent. It has been found that the ratio 
02/01 from an individual fit to the CLEO data of B — > D^n(p) || varies 
significantly ||. It has even been shown that an allowed domain (ai, 02) exists 
for the three classes of decays B° -> Z>W+, B° -> _D ( * )0 and B~ -> DW°, 
only when the experimental errors are expanded to a large extent ||. On 
the other hand, a negative 02/01 and a positive 02/01 were concluded from 
the data of charm and bottom decays ||], respectively and this subject 
remains controversial. 

In the BSW model nonspectator contributions from P^-exchange (or an- 
nihilation) diagrams are not included. In previous studies such contributions 
were assumed to be negligible, though a convincing justification is not yet 
available. Hence, it has been a challenging task to evaluate W -exchange 
contributions to nonleptonic B meson decays reliably. 

Recently, we have shown that the perturbative QCD (PQCD) formalism, 
which contains the resummation of large radiative corrections, is applicable 
to B — > D decays in the fast recoil region of the D meson ||, and have 
employed this formalism to investigate the spectator contributions from ex- 
ternal iy-emission diagrams to B -> £> w decays 0. In our approach all the 
relevant form factors can be evaluated explicitly without resort to models. 
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For neutral B meson decays such as B° — > D^ + n~ (p~), there are additional 
contributions from J^-exchange diagrams. For charged B meson decays such 
as B~ — > D^°n~ (p~), both external and internal VT-emission diagrams con- 
tribute. 

In this letter we shall extend the PQCD formalism to the study of internal 
W^-emission and VT-exchange diagrams, whose contributions are computed 
reliably. Our work helps clarify the ambiguity in the extraction of a\ and a 2 
in the BSW model. A PQCD analysis of these contributions based on the 
Brodsky-Lepage theory [Q] has been performed in || . However, this approach 
was criticized in The parametrization of the fractional momentum car- 



ried by the light valence quark in a B meson does not satisfy the on-shell 
requirement from the parton model, such that extra imaginary contributions 
are introduced. Moreover, the predictions depend strongly on the model of 
meson wave functions [§, [TOfl . 

2. Resummation 

We review the investigation of the external jy-emission diagrams as shown 
in Fig. 1(a), which involve six form factors i — +, — , V, A\, A 2 and A 3 , de- 
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fined by the hadronic matrix elements of vector and axial vector currents M. 
Assume that the momentum P 1 (P2), the mass M B (M fl (.j) and the velocity 
vx (v 2 ) of the B (DW) meson are related by Pi = M B v 1 (P 2 = M dM v 2 ). The 
velocity transfer 77 — V\ • v 2 is expressed in terms of the momentum transfer 
q 2 = (Pi - P 2 ) 2 as 

' 2M B M DM ' 1 ; 

In the infinite mass limit of M B and M D (*), the six form factors have the 

relations 

e+ = = u = u = e, e- = a = 0. (4) 



£ is the so-called Isgur-Wise (IW) function |TTJ , which is normalized to unity 
at zero recoil rj — > 1 by heavy quark symmetry (HQS). 

In the rest frame of the B meson Pi is written as Pi = (Mb/\/2)(1, 1, Oy), 
and P2 has the nonvanishing components 

f/^**, (5) 

Let fci (^2) be the momentum of the light valence quark in the B (D^*>) 
meson, satisfying k\ pa ik\ ~ 0). fei may have a large minus component 
fcf , defining the momentum fraction x\ = frf/Pf", and small transverse com- 
ponents kiy. ^2 may have a large plus component k 2 , defining x 2 = k 2 jP 2 ', 



and small k 2 <r- As r] — > 1 with P 2 + = P2 = / v2, the two meson wave 
functions strongly overlap, and the form factors £j are dominated by soft con- 
tributions. In the large rj limit with P 2 + 3> M D (»> / v^2 ^> P 2 ~, the £)W meson 
behaves like a light meson ||. Then P — > P/*) transitions occur through 
hard gluon exchanges, to which PQCD is applicable. 

Double logarithms from radiative corrections to the B (P^) meson wave 
function <pB (0d(*)) have been organized into an exponent s using the resum- 
mation technique 0. We quote the results as 



^(xi,PiA,/i) = B (a;i)exp 



-s{x 1 P{M) ~ 2 / —j(a s (p)) 

JX/bi A* 



!> D w(x 2 ,P 2 ,& 2 ,/i) = jD (*)(x 2 )exp 



dfi 

l/bi p, 

-s(x 2 P 2 + ,6 2 )-s((l-x 2 )P 2 + ,6 2 ) 



-2 f — 7(« s (Al)) 
•/1/&i> u 



(6) 



where b\ (6 2 ) is the Fourier conjugate variable of k\x (P 2 t), and can be re- 
garded as the spatial extent of the P (PW) meson, [i is a renormalization 
and factorization scale. 7 = —a s /ir is the quark anomalous dimension. The 
expression of s is very complicated and referred to [^, The initial con- 
ditions 4>i(x), % = B, D and D*, are of nonperturbative origin, satisfying 
the normalization Jq 4>i(x)dx = /j/(2\/6), being the corresponding meson 
decay constants. 
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The evolution of the hard scattering amplitude H is expressed as || 

H(ki, k£, b h b 2 , /i) = H(ki, k£, bi, b 2 , t) exp 

where the variable t denotes the largest mass scale of H. Combining Eqs. (@) 
and (|7|), we obtain the factorization formulas of £j. 

For nonleptonic decays, there exist additional important corrections from 
final-state interactions with soft gluons attaching the outgoing hadrons. It 
has been argued that these corrections produce only single logarithms |J, 
and are thus not considered here. 



J fx n 



(7) 



3. Nonleptonic B —> decays 

In the analysis of nonleptonic B meson decays we employ the Hamiltonian 
in Eq. instead of the effective Hamiltonian in ([[]). We shall address how 
to formulate our PQCD theory based on the effective Hamiltonian in the end 
of this letter. The amplitudes for charged and neutral B meson decays are 
then written as 



^=v cb v: d 



7i-\{du)\0)(D^ \{cb)\B-) - (7T'\(D wu \(du)(cb)\B~) 



(*)0| 



( ' F V cb V: d \(n~\(du)\0)(D^\(cb)\B } - (n-\(D^\(du)(cb)\B } 



V2 



(9) 



The first terms in Eqs. (|8]) and (|9]) correspond to the external W^-emission con- 
tributions, and the second terms to the internal VT-emission and W -exchange 
contributions, respectively. Note the minus signs in front of the second terms, 
which are associated with the interchange of two u (d) quarks in the inter- 
nal ^-emission (W -exchange) diagrams as shown in Fig. 1(b) (Fig. 1(c)) 
compared to Fig. 1(a). There are another two nonfactorizable internal W- 
emission diagrams, in which one end of the gluon line attaches the quark 
of the D^*' meson. However, it can be shown that the contributions from 
these two diagrams cancel partially, and are thus less important. For the 
same reason, we neglect the contributions from another two nonfactorizable 
VT-exchange diagrams, in which the gluon attaches the quark of the B meson. 

It has been argued that the internal H^-emission and W-exchange am- 
plitudes can be evaluated in the PQCD approach [|7|. The decay rates of 
B —> D(*> transitions have the expression 

= -^Gl\V cb \ 2 \V ud \ 2 ml^^\M t \ 2 , (10) 
1287T r 

with r = M DM /M B and i = 1, 2, 3 and 4 denoting B~ -> D tt", B° -> 
D + n~, B~ — > D*°ti~ and B° — > D* + n~, respectively. The amplitudes Mi 
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are written as 



Mi 
M 2 
M 3 
M 4 



f. 



M(l + r)e+ - (1 - r)e_] + ^mt , 



7T 



/ T [(l + r)( + -(l-r)(_] + -M c 



(11) 
(12) 



1 + r 

2r 
1 + r 

2r 



/,[(! + r)U - (1 - r)(rU + U)] + ^-&t , (13) 



7T 



/.[(I + r)U - (1 - r)(rU + U)l + 4^ ■ ( i4 ) 



The color-suppressing factors l/iV c have been shown explicitly. 

The factorization formulas for £j, i = +, V, Ai and A 3 (^ 3 = £y) and 
for £j, j = — and A 2 are given by 



167rCirA/rM| / dx\dx 2 / bidb 1 b 2 db 2 (j) B (xi)(j) D {*)(x 2 )a s (t) 
Jo Jo 

x [(1 + (iX 2 r)h(xi,x 2 , h, 6 2 , m) + (r + (-xi)h(x 2 , x\, 6 2 , 6i, ra)] 

x expf-S^i, ^2, &i, 6 2 )] , (15) 

IQuCfVtMq / dx\dx 2 \ bidbib 2 db 2 (j) B (xi)(fi D (*)(x 2 )a s (t) 
Jo Jo 

x [Qx 2 rh(xi, x 2 , bx, b 2 , m) + (' j xih(x 2 , Xi, b 2 , b 1: m)] 



x exp[-S , (xi,x 2 , 6i,6 2 )] 



(16) 



with the constants 



c + = c; 



.3 

2 



/77 — 1 

77 + 1 



C- = -C- 

Ca 2 = o , Cm = -i 



77 + 1 
1 

~2 



' bAl 2(77 + 1) 
1 [rj+ 1 / 3 



77 2 ' V *7 - 1 
V 



"-2 



(17) 



Cp = 4/3 is the color factor. These form factors are evaluated at the maximal 
recoil 77 = ?7 max = (1 + r 2 )/(2r) in Eqs. fln|) to ([141) . 
The form factors and are given by 



A*) 

^int 



*(*) 

S^cxc 



167rCF\/rM| / dx x dx z / M&iM&30B(zi)07r(z3)a«(*int) 

J JO 

x (1 + a? 3 (l - r 2 ))/i int (xi, x 3 , 61, 63, mint) 
+Cint )a; i r2 ^int(a;3, ari, &3, &i, ^int) exp[-Si n t(si, £3, &i, 63)] , (18) 

rl roo 

16nC F y/rMl / dx 2 dx 3 / b 2 db 2 b- i db- i (l) D {,)(x2)(l) 7 ,{x^)a s (t exc ) 
Jo Jo 

X (JC 3 (1 - r 2 ) - Cexc r2 ) /l exc(^2, S3, 6 2 , &3, ™exc) 

— ^2^exc(s3) x 2 , 63, 62, m e xc)] exp[-S , cxc (x 2 , X 3 , 6 2 , 63)] , (19) 

with the constants £ int = — Cint = 1 an d Cexc = — C xc = 1- In the derivation 
of we have assumed that k\ has a large plus component kf = X\P± . 
Here £3 is the momentum fraction associated with the pion, and 63 can be 
regarded as the spatial extent of the pion. and £W are in fact the B — > 7r 
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and D — > 7r transition form factors, respectively, evaluated at r\ = r] max . 

The Sudakov exponents S"s, which group the exponents in Eqs. (|6|) and 
(0), are given by 



5 = s(x 1 P 1 -,b 1 )+s(x 2 P+,b 2 ) + s((l-x 2 )P+,b 2 ) 



1 

"A 



In Ht / A l I In ln( ^ /A) 



5- 



int 



■ln(6iA) ' "~-ln(6 2 A) 
s( Xl P+, h) + s(x 3 P 3 , b 3 ) + s((l - z 3 )P 3 ~ 6 3 ) 



1 



ln(t int /A) ln(t int /A) 
In — ; — - — — + In 



-ln(6iA) -ln(6 3 A)_ 

s{x 2 P 2 + , b 2 ) + s((l - x 2 )P 2 + , b 2 ) + s(x 3 P 3 ", 6 3 ) 



+s((l-x 3 )Ps,b 3 



1 

ft 



ln ln(t exc /A) ^ ln(t exc /A) 



ln(6 2 A) 



ln(fo 3 A) 



(20) 



(21) 



(22) 



with /3i = (33 — 2n^)/12 and rij = 4 the number of flavors. The QCD scale 
A = Aqcd w iU be set to 0.2 GeV below. The factors e~ s fall off quickly in 
the large b, or long-distance, region, giving so-called Sudakov suppression. 
In Eqs. fD, (T 



|), (|lq) and ([L9]) the functions /i's, obtained from the 
Fourier transform of the lowest-order H, are given by 



h(x 1 ,x 2 ,bi,b 2 ,m) = K (^XiX 2 mbi 



x [6{b x - b 2 )K ( y /x 2 mb 1 ) I {y/x 2 mb 2 ) 

+6{b 2 -b 1 )K {y/x^Fib 2 )I {^ 2 !Eb 1 )] , (23) 
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h- m t(xi,x 3 ,b 1 ,b 3 ,m in t) = h(xx, x 3 , b u b 3 , m- m ) , (24) 
h exc (x 2 ,x 3j b 2j b 3l m exc ) = H§> ( y/x 2 x 3 m exc b 2 ) 

x 6(b 2 - 6 3 )#o {\/x 3 m exc b 2 ) J (^x 3 m exc b 3 ) 
+9(b 3 - b 2 )H^ ] (^x 3 m exc b 3 ) J (^/x 3 m exc b 2 ) , 

(25) 

with m = (77 + y/rj 1 — 1)M B M D (*) and mi nt = m exc = — M^ w . It is 
obvious that the W^-exchange contribution is complex due to the exchange 
of a time-like hard gluon. In the above expressions the large scales t's take 
the values 

t = max(y/xix 2 m, l/&i, I/&2) (26) 
t int = max( v ^m^, l/6 3 ) (27) 
t exc = max(y/x 2 x 3 m exc , 1/&2, 1/&3) • (28) 

We choose f B = 200 MeV, f D = f D * = 220 MeV §, = 0.043 for the 
CKM matrix element || |7j, and the Chernyak-Zhitnitsky model [|TI| 



^(*) = ^A*(l-x)(l-2:r) 2 (29) 

for the pion wave function, /„. = 132 MeV being the pion decay constant. For 
the B meson, we employ the wave function from the relativistic constituent 
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quark model [ O , 



>B{X,k T 



c B + M ° 



+ 



1 -2 



1 — x x(l — x) 



(30) 



The normalization constant Nb and the shape parameter C B are determined 



by the conditions 



dx 



dx 



d 2 kn 



167T 3 
T 



d 2 h 



1Qtt' c 



(f) B (x, k T )] 2 



B 



2VQ 
1 

~ 2 ' 



(31) 



The B meson wave function is then given by 



d 2 kn 



A 7 , 



x(l — x) 



(32) 



167r 3^~>"^ 167t 2 M| + Cb(1 -x) ' 
with N B = 650.212 and C B = -27.1051. We assume that the L> w meson 



wave function possesses the same functional form as Eq. 



(j) D {*){x) 



N DM x{\ - xf 
16vr 2 M 2 DM +C DM (l-x) 



(33) 



but the shape parameter C$ is fixed by the B — > data 



4. Discussion 



We adopt G F = 1.16639 x 10~ 5 GeV" 2 , \V ud \ = 0.974, M B = 5.28 GeV, 



M D = 1.87 GeV, M D * = 2.01 GeV 0, and r B o = 1.53 (r B - = 1-68) ps for 



13 



the B° (B ) meson lifetime [T5]. The experimental data of the branching 
ratios are B(B~ -> D°tt-) = (5.5 ± 1.1) x 10~ 3 , B(B° -> D+vr-) = (2.9 ± 
1.2) x 10- 3 , -> D*%-) = (5.2 ± 1.7) x 10~ 3 , and B(B° -> D*+7r-) = 

(2.6 ± 0.8) x 10 -3 0. We determine the parameter , ie., the D^> meson 
wave function, by fitting our predictions to B(B° — > D + n~) and i3(5° — > 
D*+7T-). They are found to be C D = -2.95 GeV 2 and C D . = -3.05 GeV 2 . 
The corresponding normalization constants are then = 128.77 GeV 3 and 
N D » = 173.48 GeV 3 . 

With these parameters we derive the branching ratios B(B~ — > D°tt~) = 
4.41 x 10~ 3 , B(B° -> D+TT-) = 2.91 x 10" 3 , B(B~ -> D*°tt-) = 3.97 x 10~ 3 , 
and S(5 — > .D* + 7r~) = 2.59 x 10 -3 , which are within the errors of the data. 
It is observed that the internal VT-emission amplitude is additive to and about 
20% of the external jy-emission amplitude. Therefore, our analysis favors 
a positive 02/01 for B meson decays. The internal W^-emission contribution 
obtained in || is complex, but it is real in our formalism, and is twice larger 
in magnitude. Hence, our predictions for the charged B meson decays are in a 
better agreement with the data. We also find that the W-exchange amplitude 
is mainly imaginary, and its magnitude is only 3% of the external jy-emission 
amplitude. The real part of the iy-exchange contribution is at least one 
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order of magnitude smaller than the imaginary part. The corresponding 
results in |J have the same order, but the real part is larger or only slightly 
smaller than the imaginary part in magnitude for different models of pion 
wave functions. The branching ratios of the charged B meson decays are then 
(1.2) 2 x (tb-/tbo) ~ 1.5 times of those of the neutral B meson decays, fairly 
consistent with the data, and with the prediction B(B~ — > D^*'°tt~)/B(B° — > 
DW+tt") « 1.7 from the BSW method @. 

The branching ratios of the charged B meson decays fall below the central 
values of the data. The agreement can be improved by including the two 
nonfactorizable internal VK-emission diagrams. Because of the complexity of 
their evaluation, we shall not study these two diagrams in this letter, but 
discuss them elsewhere. 

It is worthwhile to exhibit the spectrum dT / dq 2 of the semileptonic decay 
B° — > D* + £~v. For this process only the external ^-emission diagrams 
contribute, and the expression is written as [0] 



x ^ 2(1 - 2rjr + r 2 ) 



e A M + v —\ev(v) 
77+1 



+ [(v - r)UM -(v-i) HaM + UM)] 2 \ , (34) 
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where the form factors £i(r}) have been defined in Eqs. (|I5D and (|16l). It is 
observed from Fig. 2 that our predictions match the data |T5| at low q 2 , but 
begin to deviate above q 2 = 4 GeV 2 , the slow recoil region in which PQCD 
is not reliable. 

Since the W-exchange contributions are negligible, and the B — > 7r and 
B — > p form factors are roughly equal [|16|], the branching ratios of the de- 
cays B — > D^p can be easily estimated by substituting the p meson decay 
constant f p = 220 MeV for /„. in the corresponding formulas. Assuming a 
vanishing p meson mass, we have B(B~ — ► D°p~) = 1.2%, B(B° — > D + p~) = 
8.1 x 10" 3 , B(B~ -> = 1.1%, and -> = 7.2 x 10" 3 . 

Compared to the data (1.35 ± 0.30)%, (8.1 ±3.6) x 10" 3 , (1.68 ±0.58)%, and 
(7.4 ±2.7) x 10~ 3 [U, respectively. Our predictions are satisfactory except for 
the mode B — > D*°p~. Similarly, the overall consistency can be improved 
by including the two nonfactorizable internal jy-emission diagrams. 

It is nontrivial to incorporate the effective Hamiltonian in Eq. (|J) into 
factorization theorems. The Wislson coefficients c\ and C2, or a\ and a,2 
equivalently, which take into account the evolution from the W boson mass 
My/ down to the scale /x, arise from the irreducible radiative corrections to 
the decay vertex. In the conventional approach \x is set to a value of order 
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Mb, and thus a scale dependence is introduced. We argue that /i in a\ and 
a 2 should be set to the large scale t of the hard scattering. The scale t is 
then integrated over, such that predictions remain RG invariant. Hence, a\ 
and CZ2 are not parameters as in the BSW model, but vary according to their 
evolutions in our approach. The controversy over the scale setting in the 
analysis of inclusive nonleptonic B meson decays may be resolved by our 
RG invariant formalism. 

An immediate observation from the above choice of fi is that the internal 
VT-emission amplitudes are still positive in bottom decays. If main con- 
tributions came from the region with t close to the b quark mass Mb, the 
coefficients of the external and internal jy-emission amplitudes, a\ = 1 and 
a 2 = l/N c , are replaced by ai(M b ) = 1.03 and a 2 (M b ) = 0.11. The inter- 
nal jy-emission amplitudes, however, may become negative for charm decays 
because of ai(M c ) = 1.09 and a2(M c ) = —0.09, M c being the c quark mass. 
Certainly, this issue on the sign of a 2 /cti [Jj, |B| needs more quantitative study. 
The details will be published in a separate work. 

Effects from nonfactorizable final-state interactions can be included 
into our formalism easily, which will be discussed elsewhere. The discrepancy 
between model-dependent theoretical predictions and experimental data as- 
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sociated with the decays B — > ipK^ |T8[ can also be investigated. 
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Figure Captions 

Fig. 1. (a) External W^-emission, (b) internal VT-emission, and (c) W- 
exchange diagrams with the b quarks and the W bosons represented by double 
lines and dashed lines, respectively. 

Fig. 2. The spectrum dT/dq 2 of the semileptonic decay -B — > D* + £~u. 
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(c) 
Figure 1 




Figure 2 



